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One of the key steps in nanotechnology is our ability to engineer and fabricate low-dimen- 
sional nano-objects, such as quantum dots, nanowires, two-dimensional atomic layers or 
three-dimensional nano-porous systems. Here we report evidence of nanotunnel opening 
within the subsurface region of a wide band-gap semiconductor, silicon carbide. Such an 
effect is induced by selective hydrogen/deuterium interaction at the surface, which possesses 
intrinsic compressive stress. This finding is established with a combination of ab-initio 
computations, vibrational spectroscopy and synchrotron-radiation-based photoemission. 
Hydrogen/deuterium-induced puckering of the subsurface Si atoms marks the critical step in 
this nanotunnel opening. Depending on hydrogen/deuterium coverages, the nanotunnels are 
either metallic or semiconducting. Dangling bonds generated inside the nanotunnel offer a 
promising template to capture atoms or molecules. These features open nano-tailoring 
capabilities towards advanced applications in electronics, chemistry, storage, sensors or 
biotechnology. Understanding and controlling such a mechanism open routes towards 
surface/interface functionalization. 
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Silicon carbide (SiC) offers fascinating structural, thermo- 
mechanical, electronic and chemical properties with a vast 
range of advanced applications including high-power, high- 
frequency and high-temperature electronic devices and sensors^"^. 
In addition, SiC has a remarkable biocompatibility making it 
useful for biomedical applications^. Furthermore, it is a substrate 
especially suitable for the growth of epitaxial graphene with very 
promising potential applications in electronics and spintronics . 
SiC exists in > 170 polytypes including cubic (3C or P), hexagonal 
and rhombohedral forms . Together with hydrogen and oxygen. 
Si and C, including SiC, are among the most abundant species in 
the universe, primarily (80%) composed of the 3C-SiC polytype. 
Indeed, at the end of the nineteenth century, Henri Moissan found 
SiC on a meteorite in the Diablo Canyon, Arizona^. The cubic 3C- 
SiC(OOl) surface is of special interest with its atomic structure and 
many of its properties driven by surface strain/ stress relief^"^^ 
This leads to a large number of surprising and interesting surface 
arrangements with more than ten different well- characterized 
surface reconstructions (such as 3x2, 8x2, 5x2, 7x2, 
(2n+ 1) X 2; c(4 x 2), c(2 x 2)), ranging from Si-rich to C-rich 
systems including notably in the latter case, graphene^"^^. The 
atomic structures of these reconstructions are now well understood 
and established on the ground of advanced experimental 
techniques such as scanning tunneling microscopy (STM), 
synchrotron- radiation-based photoemission spectroscopy/ 
diffraction (SR-PES, SR-PED), grazing incidence X-ray 
diffraction, low-energy electron diffraction, near-edge absorption 
fine structure and state-of-the-art ab-initio calculations^"^^. It is 
therefore important to deepen our understanding of this material 
so that its full potential can be realized. For these reasons, SiC and 
its surfaces have attracted special attention and a range of exciting 
and unique properties have been discovered. Among many others, 
the self-formation of highly stable massively parallel passive or 
active atomic lines and nanowires at the surface^' and the 
interaction of H/D atoms with the Si-rich 3C-SiC(001)-3 x 2 
surface reconstruction are of special interest^^"^^, with the 
first example of H/D -induced metallization of a semiconductor 
surface^^'^^. 

The atomic structure of the Si-rich 3C-SiC(001)-3 x 2 surface 
reconstruction is now very well understood and established on 
the basis of atom-resolved real-space STM^, synchrotron 
radiation-based grazing incidence X-ray diffraction^ and PED^^ 
experiments, and ab-initio calculations . The structure includes 
three Si atomic layers at 1/3, 2/3 and 1 ML, respectively, on top of 
the first C atomic plane^"^^ as illustrated in Fig. 1. The first plane 
(1/3ML) includes dimer rows having asymmetric Si-Si dimers all 
tilted along the same direction, the second plane (2/3 ML) also 
includes dimer rows having alternating lengths and the third 
plane also has dimer rows^~^^. This complex structure tends to 
reduce the surface strain. Notice that the H(D)-induced 
metallization of the 3C-SiC(001)-3 x 2 surface takes place while 
its initial 3x2 array is maintained^^'^^. 

The present investigation provides evidence of 'nanotunnel' 
formation. Unlike atomic lines or nanowires that are self- grown 
on a surface such a nanostructure is not generated by self- 
assembling of atoms, but by the auto-organization of nano-voids 
below the surface. This void structure is induced by the 
adsorption of atomic hydrogen/deuterium (H/D) on a 3C- 
SiC(001)-3 X 2 reconstructed surface. The critical effect, which 
opens such nanotunnels, occurs by H- induced puckering of Si 
atoms in the subsurface layers connected with stress relief, 
thereby decoupling the second Si layer and activating the third Si 
layer below the surface. These nano-scale void structures are well- 
oriented and periodically spaced, offering unique opportunities 
for surface functionalization. In the literature, there are examples 
of voids generated at a surface or below. However, they are in the 



jim or sub-jim or mesoscopic scales and exhibit no ordering or 
self- organization mechanisms, since they are generated by 
experimental techniques such as ion bombardment or 
etching^^"^^. Nano-cavities generated by ion bombardment of a 
copper surface are, so far, the only voids known at the nanometre 
scales, but again, they exhibit no spatial ordering^"*. 

In this work, we use high-resolution electron energy loss 
spectroscopy (HREELS) and synchrotron-radiation-based photo- 
emission spectroscopy (SR-PES) combined with ab-initio calcula- 
tions to reveal the opening of nanotunnels induced by H(D) 
exposures of a 3C-SiC(001)-3 x 2 surface. Nanotunnel opening 
activates the Si atoms in the third layer inside the tunnel. H atoms 
binding to these Si atoms transform the system from a 
semiconducting to a metalHc state. Further exposure returns the 
system to a semiconducting status. The nanotunnel structure is 
consistent with results from SR-PES, HREELS and earlier infrared 
absorption spectroscopy (IRAS), STM and SR-based core-level 
PES results^^'^^'^o. 

Results 

Atomic structure and electronic properties by ab-initio 
simulations. We first look at the results of ab-initio calculations 
for the 3C-SiC(001)-3 x 2 hydrogenation as illustrated in Figs 1 
and 2. Starting with the reconstructed 3x2 clean surface, H(D) 
atoms first bind to the Si dangling bonds in the first layer, 
denoted by Si la and Si la' in Fig. 1 and seen by STM^^ The 
bonding is strongly exothermic releasing 3.23 eV per H atom. The 
resulting structure, labelled as '2 H atoms' is characterized by a 
symmetric surface dimer with the dangling bonds saturated by 
H(D) atoms, in excellent agreement with STM results^^. In the 
next step, H atoms are placed in the grooves between the dimers 
and the atom positions are relaxed by energy minimization. The 
calculations show that the Si atoms in the second layer, labelled 
Si2a, trap H atoms with an energy release of 1.50 eV per atom. 
This step breaks the bonds between Si atoms in the second and 
third layer (Si2a and Si3a), leading to a puckering around the Si2a 
atoms, and the opening of a nanotunnel. This rearrangement 
produces reactive Si atoms in the third layer (Si3a). The 
unsaturated bonds of the Si3a atoms give rise to a double peak 
in the density of states (DOS) between 0 and — 1 eV as shown for 
the 6 H structure in Fig. 2. The surface remains semiconducting. 
The formation of the nanotunnel opens the access to the third- 
layer Si3a atoms while sterically hindering the path to the Si3b 
atoms. Hence, additional H atoms are likely to bind to the empty 
dangling bond of the Si3a atoms in the third layer, forming a 
bonding arrangement H-Si-C-, since the third-layer Si atoms are 
bonded to C atoms. A partial saturation of the Si3a atoms with H 
atoms causes a partially filled band at the Fermi level (see blue 
parts of the DOS in Fig. 2). In other words, at this H-coverage the 
surface appears to be metallic. Additional exposure to H saturates 
all dangling bonds of the Si3a atoms and the surface returns to a 
semiconducting state. 

As will be discussed in the following, the computed vibrational 
frequencies for these configurations as well as the semiconducting- 
metallic-semiconducting sequence upon increasing H exposure are 
perfectly consistent with our experimental results for vibrational 
and electronic properties. 

Vibrational properties measured by HREELS. Next, we look at 
the HREELS measurements. Compared with IRAS, HREELS 
extends over a wider frequency range and can be performed at 
low sample temperatures, offering the opportunity to explore the 
low frequency modes predicted by theory. Figure 3a,b display the 
elastic peaks for the clean and SOL (1L=1 Langmuir = 10~^ 
torrs) H-exposed 3C-SiC(001)-3 x 2 surfaces at 300 and 105 K, 
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Figure 1 | Structures of the clean and H/3C-SiC(100)-3 x 2 surfaces and views of a nanotunnel. (a) The 2 H, 6 H, 8 H and 10 H structures refer 
to the number of H atonns per surface dimer. Nanotunnel opening is seen for 6 H, 8 H and 10 H. The energy release for each step is indicated, as, for 
example, 2( — 3.23) eV, that is, 3.23 eV per H atom. Si1a(b), Si2a and Si3a(b) refer to Si atoms in the first, second and third plane, respectively. Si dangling 
bonds are depicted in blue, (b) The nanotunnel has a size of 6.2 x 4.4 A (at the level of the second layer). The red atoms within the nanotunnel 
mark Si atoms with one or two dangling bonds corresponding to the 8 H (metallic) and 6 H (insulating) structures, respectively. 



respectively. The elastic peak full-width-at-half-maximum of the 
H-exposed surface is broadened^^. This broadening is 
significantly larger at 300 K than at 105 K, indicating either 
unresolved multiple low frequency plasmon excitation or the 
onset of a Drude damping^^ (see also the HREELS experiments 
section in Methods). In both cases, this means that the 3C- 
SiC(001)-3 X 2 surface exposed to SOL H has a metallic character. 

We now compare the vibrational frequencies measured by 
HREELS for the SOL H(D)/3C-SiC(100)-3 x 2 surface to those 
computed for the 8 H metallic model that are displayed in Table 1 
together with previous IRAS^^ measurements. We also investigated 
the deuterated system to explore regions hidden by the Fuchs- 
Kliewer (FK) modes. The corresponding vibrational spectra for 
H(D) are displayed in Fig. 4a,b for in and out of specular 



conditions, respectively, providing the description of the modes- 
see Table 1 also listing previous IRAS results^ . Measurements were 
performed at 300, lOS and 30 K obtaining very similar results. The 
losses at 940, 1,880 and 2,820 cm ~^ correspond to single and 
multiple excitations of the SiC FK mode^^, and are also present on 
the bare surface. The smaller losses at 63Scm~^ moving at 
600 cm ~ ^ out of specular (43Scm~ ^ out of specular only) and at 
2,130 cm ~^ moving to 2,120 cm ~^ (1,S40-1,S32 cm~ ^ for the 
deuterated surface), are due to the bending and stretching of the 
H(D)-Si bonds, respectively. The latter region had been previously 
explored also with IRAS at room temperature only, finding 
compatible results^^. 

These spectral features are found to be in remarkable 
agreement with theory accounting for all observable modes 
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Figure 2 | Atom-projected DOS for nanotunnels. The atom-projected DOS are displayed for the nanostructures corresponding to the (a) 6 H, (b) 8 H and 

(c) 10 H nanotunnels shown in Fig. la. For the metallic 8 H nanotunnel (b), we notice the DOS build-up at the Fermi Level £'f (shown in blue), which is 
occuring predominantly in the third Si layer for the Si3a atoms. Note also the semiconducting-metallic-semiconducting transition upon increasing H 
coverage, with H atom-induced partially empty Si dangling bonds (shown in blue in Fig. 1) for the metallic structure. 
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Figure 3 | HREELS specular elastic peak for clean and SOL H exposed 3C- 
SiC(001)-3 X 2. (a) HREELS spectra recorded at a temperature of BOOK 
and (b) at a temperature of 105 K. The factor 2 broadening observed at 
300 K for the hydrogenated surface is indicative of surface metallization. 
Such broadening is not present at 105 K in agreement with the suppression 
of the mechanisms of either multiple excitation of a low frequency plasmon 
or Drude damping. 

involving H(D) vibrations. The measured frequency of 
635(435) cm ~ ^ is associated with a symmetric wagging mode 
of the Si2a-H(D) group, possibly with an admixture of the 
wagging mode of the Sila-H(D) atoms. Owing to the inclination 
of the Si2a-H(D) bond towards the surface plane, the wagging 
mode has a strong component perpendicular to the surface and 
thus is observable by HREELS, while the wagging of the Sila- 
H(D) is predominantly parallel to the surface resulting in a much 
weaker signal. The computed frequencies of 600 (438) cm ~ ^ for 
Si2a-H(D) and 621 (422) cm for the Sila-H(D) modes agree 
with experiments within the accuracy of vibrational modes 
obtained by density functional theory (DFT) -generalized gradient 
approximation and the harmonic approximation. In fact, 
anharmonicities can change frequencies up to ~ 80 cm ~ ^ for 
groups such as surface silanols^^, which provides an indication of 
the accuracy of the present calculations. As evident from Table 1, 
the largest difference between experimental and theoretical values 
is 37cm~^ with the theoretical value always lower than the 
experimental one, indicating the systematic nature of this error. 



Two high frequency modes have been observed with IRAS at 
2,140 and 2,1 18 cm ~^ (ref. 23). The upper frequency mode has 
been associated with a Si-H stretch for Si bonded to C, whereas 
the lower mode is related to H on the topmost Si dangling 
bond^^. The former has a higher frequency due to charge transfer 
from Si to C, which slightly enhances the ionic character of the 
Si(^^)-H(^~) bond compared with a Si-Si-H bonding 
arrangement as it is the case for the Sila-H bond, in agreement 
with the higher electron affinity of C compared with Si. The 
separation of the two modes by 22 cm ~ ^ is very well captured by 
the present calculations for the 8 H structure, which show modes 
at 2,120 and 2,093 cm ~^ whereas HREELS gives 2,130 
(1,540) cm ~^ for the symmetric stretch of the surface dimer, 
which may be influenced by signals from the third layer. Off- 
specular HREELS data may provide insights on the anti- 
symmetric stretch modes too, and show indeed slightly lower 
frequencies. In agreement with experiment, the computations 
give a lower frequency for the anti- symmetric modes compared 
with the symmetric stretch ones of the Sila-H bond. The 
computed Si2a-H(D) stretch mode at 2,020 (1,457) cm ~ ^ is not 
observed owing to its weak out-of-plane component. Thus, all 
observed H(D) experimental modes are accounted for by the 
vibrational properties of 8 H structure. One should note that the 
modes associated with atoms Sila and Si2a of the 10 H structure 
are similar to those of 8H. No modes are present at the 
frequencies of 778 cm ~^ for D-exposed or 1,450 cm ~^ for 
H-exposed 3C-SiC surface as predicted by some of the DFT 
calculations for vibrations involving H(D) in bridge positions in 
the third plane^^"^^. Modes at 1,100 cm ~^ (H) and 1,025 cm ~^ 
(D) are undetectable due to the overlap with the tail of the strong 
FK mode at 940 cm ~^ In particular, we note that the 
corresponding vertically polarized mode at 1,450 cm ~^ would 
be easy to observe since it would nicely faU in-between the FK 
peaks, which is not the case. This excludes H(D) atoms in a bridge 
bond configuration as a suitable model to explain the metaUicity as 
suggested in some DFT calculations^ Extended to low crystal 
temperatures, our HREELS investigation demonstrates no 
subsequent structural change. Also, the H interaction with the 
second layer has been invoked in refs 29,31,32. A recent 
calculation suggested that H(D) atoms binding to the 2nd Si 
layer rather than to the 3rd one in bridge positions can explain the 
metallization^^. However, this study falls short in explaining the 
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Table 1 | Measured and computed vibrational frequencies of the metallic structure. 
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The vibrational frequencies measured (SOL H(D)-exposed) and computed (8 H(D)) for the metallic structure (see Fig. 1) are shown above. Frequencies are in cm ~ I The modes are strongly coupled. The 
HREELS data (SOL H/D) recorded out of specular are sensitive also to the antiphase motion of the H(D) atoms bonded to first and third Si layer, while only the symmetric modes can be excited in- 
specular The computed values in parentheses for the Si(la)-H(D) stretch mode are for a structure with H(D) atoms bonding only to Sila atoms of the surface dimer. 
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Figure 4 | HREELS and computed vibrational frequencies for H/D-induced metallic nanotunnel (a) HREELS spectra recorded at a temperature 7=105 K 
for H (blue dots) and D (red continuous line) in specular and (b) out of specular conditions. The spectrum for the clean surface (black dash) is 
shown for comparison, too. Note the combination band between the FK and H(D) stretch at 3040 cm ~^ (2,460 cm ~^). The spectra correspond 
to SOL H(D) exposed 3C-SiC(001)-3 x 2 surfaces. The vibrational frequencies measured by HREELS for H and D are marked by ^ and ^ bars, and 
labelled in blue and red respectively, whereas the corresponding calculated frequencies are labelled in green and marked by a — — bar. Vibrational 
frequencies for the clean 3C-SiC(001)-3 x 2 surface are labeled in black and marked by a ^ bar. 



experimental data also including previous IRAS^^'^"* and core-level 
and valence-band SR-PES^'^^'^ '^^ experiments. 



Synchrotron radiation-based photoemission measurements of 
the DOS. We now focus on the corresponding electronic 
properties as measured by SR-PES and computed DOS for 
the clean and H-exposed 3C-SiC(100)-3 x 2 surface, see Figs 5 
and 6, respectively. The SR-PES spectrum from the bare 
surface^^ shows two dominant features labelled 'B' (related to 



bulk states) and 'S' (surface state related to the Si dangling 
bond)— see Fig. 5. The computed DOS (Fig. 6a) is 
integrated over the top three Si layers and the first C layer 
for the 2 H, 6 H, 8 H and 10 H structures. The photoemission 
spectra recorded at photon energy of /zv = 60eV (Fig. 6b), 
in the Fermi level (Ep) region for 2L H exposure 
(semiconducting), 36L H exposure (metallic) and 84L H 
exposure (semiconducting) surfaces are in excellent agreement 
with the calculated DOS for the 2 H, 8 H and 10 H structures, 
respectively (Fig. 6a). 
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Figure 5 | Comparison between computed and measured DOS for clean 
3C-SiC(001)-3 X 2. The calculated DOS for the clean surface is compared, 
with the experimental valence band photoemission spectrum of 3C- 
SiC(001)-3 X 2 taken from ref. 23 at a photon energy of /iv = 16.85 eV. B 
and S stay for bulk and surface contributions, respectively. Note the 
excellent agreement between theory and experiment, in particular for both 
the shape and energy position of the bulk (B) and surface (S) electronic 
states. 
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Figure 6 | SR-PES and computed DOS for H-covered 3C-SiC(001)-3 x 2 
surfaces, (a) Computed DOS for the 2 H, 6 H, 8 H atom and 10 H atom 
structures integrated over the top three Si layers and the first C layer - See 
Fig. 2. (b) SR-PES spectra recorded at a photon energy of /iv = 60eV in the 
Fermi level region for 2L H exposure (semiconducting), 36L H exposure 
(metallic) and 84L H exposure (semiconducting) surfaces. Note the 
metallic-semiconducting transition upon higher H exposures in agreement 
with the calculated densities of states. Arrows denote the position of the 
surface state for the clean SiC surface, which is quenched upon H 
adsorption. Here also, we can note the remarkable agreement between 
experimentally measured DOS and computed ones. 
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The DOS for the clean surface is close to zero at and the 
system is non-metallic as expected for a wide band-gap 
semiconductor as can be seen in the corresponding experimental 
and computed DOS in Fig. 6 (note the excellent agreement 
between the computed and experimentally measured DOS). Next, 
we now look in Fig. 6 at the effect of H adsorption. As can be seen 
in Fig. 6a,b, H eliminates already the surface state peak 'S' of the 
clean surface in the valence band at a very low 2L H exposure as 
also found in the computed DOS, with the corresponding 
structure 2 H remaining semiconducting having H atoms bonded 
to the first (Si la) and second (Si2b) atomic layers. Upon further 
H deposition, H atoms bond to the third atomic layer (Si3a) 
leading to a 8 H structure exhibiting a half-filled peak at Ep. The 
latter originates from the Si3a atom dangling bonds (see the 
computed DOS for the 8 H structure in Fig. 2) and corresponds 
to a metallic behaviour as shown in Fig. 6. Indeed, the 
photoemission spectrum shown in Fig. 6b, corresponding to a 
36L H-exposed SiC surface, exhibits a sharp spectral feature with 
a clear Fermi level indicating the metallic character of this 
structure, in excellent agreement with the calculated DOS for 8 H 
structure. Saturation of the Si3a dangling bonds eliminates these 
mid- gap states, leading the system back to a semiconducting 
status as seen for the 10 H atom structure corresponding to the 
84L H exposure (Figs 1, 2 and 6). Thus calculations and 
experiments are in perfect agreement. 

Since the synchrotron beam diameter is 300 jim with 
photoelectrons collected from an entire area covering a surface 
seven orders of magnitude larger than probed by STM/scanning 
tunneling spectroscopy (150 x 300 A)^^, evidencing H-induced 
surface metallization by SR-PES indicated that the nanotunnel 
metallization extends to a large surface. The strong metallic 
character shown in the 36L H spectrum through intense DOS and 
sharp Fermi level is clearly inconsistent with short and scarce 
nanotunnels, but in agreement with nanotunnels covering the 
whole surface (Fig. 6). 

Discussion 

The above results show a feature of central importance, namely 
the first evidence of self- organized opening of nanotunnels in a 
material. This is initiated by atomic H(D) bonding to the Si2a 
atoms, causing an inversion of the tetrahedral arrangement 
(puckering), thus breaking Si2a-Si3b bonds leading to highly 
reactive Si3a atoms. The opening of the nanotunnel frees the 
access to the Si3a atoms and additional H(D) atoms can form 
Si3a-H(D) bonds in a strongly exothermic reaction. 

At this stage, if the atmosphere contains other reactive atoms 
or chemical precursors, one can functionalize this subsurface 
layer exploiting the nanotunnel. In the present work, the surface 
is exposed only to a simple element, hydrogen, eventually leading 
to a fully hydrogenated structure. In the intermediate step, 
illustrated as a 3D view in Fig. 7, this leaves reactive dangling 
bonds inside the nanotunnel, which could be used to trap other 
species of interest. 

After nanotunnel opening, the computed properties of the 
three-nanotunnel structures are in remarkable agreement with a 
range of experimental data and there is no experimental evidence, 
which would be in contradiction to the properties computed for 
these unprecedented void nanostructures. Specifically, the points 
are as follows: (i) Agreement with HREELS: The computed 
vibrational frequencies for the nanotunnel 8 H structure with 
both H and D atoms are in very good agreement with HREELS 
data as given in Table 1 and illustrated in Fig. 4. The anti- 
symmetric stretch mode of the surface dimers as well as the out- 
of-plane wag mode of the Si2a-H(D) atoms are very well 
described. The computed Si2a-H(D) stretch mode is not 
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Figure 7 | A computed view of a 3D nanotunnei. The nanotunnel opening 
induced by the interaction of H-atoms with the 3C-SiC(100)-3 x 2 surface 
is represented for the 8 H metallic structure. Note that, in the nanotunnel, 
the dangling bonds not terminated by a hydrogen atom are not shown. 

observed, which is reasonable because this motion is mostly in 
plane and thus expected to be weak in HREELS experiments. All 
experimentally observed frequencies associated with H(D) are 
accounted for; (ii) Agreement with IRAS: The computed Si3a- 
H(D) stretch modes are higher in frequency by 27 cm ~ ^ 
compared with 22 cm ~ ^ measured in infrared spectroscopy. 
As suggested earlier^^, the stretch frequency of a Si-H bond 
should be higher if the Si atom is bonded to a C atom rather 
than another Si atom. The computed frequencies for the 8 H 
nanostructure correctly capture this situation; (iii) Agreement 
with SR-PES: The present calculations show that H atoms 
eliminate the peak originating from the asymmetric dimer on 
the clean 3C-SiC(100)-3 x 2 surface in agreement with the 
experiments. The DOS of 8 H nanostructure reveals a half- 
occupied narrow band at the Fermi level, which is associated 
with the observed metallic behaviour of this surface. Further 
exposure to H would fill this narrow band and the system 
returns to a non- metallic state, in agreement with the 
experimental observations. Furthermore, the above picture of 
H(D) atoms interacting with the third Si layer bonded to the first 
C plane is also consistent with previous core level SR-PES 
experiments indicating selective H(D) atoms interaction with 
the third Si layer^^'^^. 

These nanotunnels open a direct path for the H(D) atoms 
towards the third Si atomic plane below the topmost surface 
dimers, while impeding the access to the Si dimers located in the 
third layer as initially assumed^ ^. The ability to engineer 
semiconducting and/or metallic nanotunnels with such a novel 
one-dimensional void template offers tracks for selective 
interaction between dangling bonds generated inside the 
nanotunnel and inorganic or organic molecules, as well as 
metal atoms. After the puckering of the subsurface atoms, the 
'hole' above the activated third-layer Si atoms is relatively large, 
namely 6.2 A long and over 4.4 A wide (Fig. 1). Considering the 
covalent radius of 1 . 11 A of silicon atoms, the opening provides 
access to the nanotunnel for a large variety of atoms or molecules, 
since H is the smallest atom with atomic and covalent radii at 
0.53 and 0.31 A, respectively, leaving enough room for other 
atoms. Indeed, looking at the periodic table, we can see many 
atomic radii are smaller than 2 A (Mendeleev periodic table, 
WebElements, http://www.webelements.com/). In these views Cs, 
one of the largest atoms used, for example, in insertion 
compounds as graphite"*^ is an interesting test case. So, for 
many atoms, there is little steric hindrance for insertion into the 
nanotunnel. Regarding molecules, there are many organic and 
inorganic molecules and radicals that have size ranging below or 
near 4.4 A, small enough to make them also suitable for insertion/ 



storage within a nanotunnel - among others one can mention 
OH, O2, N2, CO, CO2, H2O, C2H2, C2H4, ethylene oxide, CH4, 
NH, NH2, NH3 and many more, which may be able to pass the 
window above a nanotunnel (Molecular Sieves, http:// 
chem.rochester.edu/~nvd/molecularsieves.html.). The critical 
geometric aspect is the 3.04 A distance between atoms Si2a and 
Si3a. The chemistry of insertion needs to be chosen such that 
there is a strong bond between the atom Si3a and the inserted 
atoms or molecules while pushing atom Si2a further away from 
the nanotunnel. 

Actually, it has already been demonstrated that a pre- oxidized 
3C-SiC(100)-3 X 2 surface could be metallized by hydrogen 
atoms, and vice-versa, that the metallization by H atoms of a 
3C-SiC(100)-3 X 2 is not removed by oxygen^'^^"^^. This leads at 
the same time to a surface with two opposite functionalization, 
metallization and passivation, a feature of very strong interest for 
interfacing with biology. Since the only remaining available 
reactive sites at the hydrogenated 3C-SiC(100) surface are the Si 
dangling bonds located within the 6 H and 8 H nanotunnels (all 
the other ones in the first and second layers being decorated by H 
atoms), this clearly indicates that the oxygen molecule/atom 
could go only into the nanotunnel. In addition, infrared 
absorption vibrational spectroscopy (IRAS) measurements have 
shown that the H-Si (2,118 cm -^) and H-Si-C (2,140 cm- l) 
bonds are not affected by oxygen exposures of the H-3C- 
SiC(100)-3 X 2 surface^"*. Therefore, one can safely anticipate that 
the same situation will occur with other molecules. 

Another interesting test case would be a metal atom like Cs, a 
'huge' atom which, despite its size, is very well known to 
intercalate easily below a surface like graphite, passing through 
the first atomic layer leading to alkali metal- intercalated graphite 
compounds^^ Indeed, the atomic and covalent radii of a Cs atom 
are 3 and 2.44 A, respectively (Mendeleev periodic table, 
WebElements, http://www.webelements.com/), making them 
very suitable to insert within a >4.4A nanotunnel having 
empty Si dangling bonds, which once again, are the only available 
'reactive' adsorption sites. It is very well known that a polarized 
covalent bonding through hybridization between the alkali 's' 
spherical orbitals with localized 'p' electronic states is taking place 
- as shown, for example, for Si^^. In addition, one should 
mention that Cs and other alkali metals are excellent catalysts, 
like, for example, of oxidation and nitridation, opening very 
interesting prospects towards selective surface functionalization"*^. 
Finally, alkali metals are known to lower significantly the surface 
work function, which could also be a very useful feature too. 

Placing metal atoms within the nanotunnels, one could also 
tune the Schottky barrier height to engineer the Fermi level (Ep) 
pinning position. Alkali metals have indeed been shown to 
selectively pin the Fermi level for other compound semiconduc- 
tors such as III-V materials like InAs, GaAs, InP or GaSb"*^""*^, 
with £p coming thereby at the highest energy level position above 
a semiconductor conduction band minimum as, for example, for 
Cs/InAs"*^. Adsorption of alkali metals have also other interesting 
properties such as their ability to control and engineer the line-up 
at heteroj unctions and insulator/semiconductor interfaces'*^. 

Other metals such as Ga or Pb could also be used as also 
shown, for example, for GaAs^^. Furthermore, inserting full 
d-band or jellium-like metals within the nanotunnel, it would 
potentially be possible to set the Schottky barrier height upon 
request, depending on the selected species as shown for, for 
example, with Au or Al on a II-VI semiconductor^^ - Note that 
the atomic and covalent radii for Ga are: 1.36 and 1.22 A, for Pb: 
1.54 and 1.46 A, for Au: 1.74 and 1.34 A and for Al: 1.18 and 
1.21 A, respectively (Mendeleev periodic table, WebElements, 
http://www.webelements.com/.), that is, once again, well below 
the 6.2 X 4.4 A hole's aperture above the nanotunnel (Fig. lb). 
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It is interesting to note that porous molecular organic nano- 
cages have been built as 'host molecules' that exhibit a selective 
interaction with incoming 'guest molecules'^^. Also, supra- 
molecular cages assembled with H atoms, have been shown to 
encapsulate molecules or metal complexes^ Similarly, 
molecules could be let into the nanotunnels either from 
the surface through the 6.2 x 4.4 A hole's aperture above the 
nanotunnel, or alternately, at a step edge, where the nanotunnels 
begin or end. The entrance could be found by shooting them 
directly in exploiting the directionality of a supersonic molecular 
beam in which the molecules are carried by an inert gas. 
As recently shown, symmetric tops can be aligned with respect 
to the rotational axis, increasing significantly the probability 
of entering with the molecular axis oriented upon purpose, 
either along or from the top, perpendicularly to the nanotunnel 

• 53-S5 

axis . 

We can thus envisage using such nanotunnels as nanotraps 
with intercalation of guest atom/molecule. Indeed, intercalation 
of atoms as large as Cs and even of Ceo molecules (having a 7 A 
diameter) under a graphene sheet has been reported^^'^^. 
Last but not least, the trapping of guest atoms can be used to 
tailor the electronic properties of the nanotunnels in analogy 
to what happens for SWCNT (single-walled carbon nanotubes) 
for which the size of the intercalated atom has been 
demonstrated to play a key role^^. The SWCNT small inner 
diameter (1.75 A) allows iodine ion insertion leading to selective 
doping^^. Once again, this indicates that a similar insertion 
should not be a problem with a nanotunnel having a 4.4 A 
width and a 6.2 A top access and 1 (8 H) or 2 (6 H) dangling 
bonds (Fig. 1). 

Therefore, insertion into the nanotunnels of most of these 
atoms or molecules could even take place through intercalation 
below a graphene layer grown or transferred onto a 3C- 
SiC(001)-3 X 2 nanotunnel surface, offering very interesting 
possibilities and opportunities to monitor the graphene/SiC 
interface including doping and its resulting properties. Indeed, a 
graphene sheet could be transferred on top of these nanotun- 
nels^^ enabling atomic scale heterojunctions or Schottky barriers 
formation. 

Interestingly, graphene layer growth can also be achieved 
using a microwave plasma technique at low substrate tempera- 
tures (300 °C instead of > 1200 °C using other processes)^^. This 
offers another promising alternative, opening engineering 
capabilities for graphene/nanotunnel interface formation 
without hydrogen atoms removal from the SiC surface/ 
subsurface or nanotunnel damage. Indeed, not only metal 
atoms like alkalis and other metals could go below the graphene 
layer (as in graphite surface), but also inorganic or organic 
molecules too, and interact with the nanotunnels as shown for, 
for example, for carbon nanotubes^^. 

There is a large variety of new physics and chemistry 
that could be used to functionalize such a SiC subsurface 
modified by hydrogen atoms, triggering interesting research 
directions, in particular in the field of graphene substrate's 
functionalization. It impacts important issues such as 
tunneling nano-junctions that are of interest in electronics/ 
spintronics applications. Using the nanotunnels for 'subsurface 
wiring' with vertical molecular connections to the outside 
could potentially be possible and presents another perspective 
that could be explored. Furthermore, compounds other than 
SiC showing the combination of surface and chemically 
induced puckering of subsurface atoms could lead to 
systems with larger nanotunnel dimensions. So, these findings 
could be useful in opening-up a class of promising 
nano-materials having a broad field of potentially interesting 
engineering capabilities. 
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Figure 8 | HREELS specular elastic peak. The elastic peak is measured for 
different crystal temperatures (30, 105 and 300 K) for (a) clean and (b) 
hydrogen-covered 3C-SiC(001)-3 x 2 surface under identical experimental 
conditions and spectrometer settings. It shows the subsequent peak 
broadening upon a 50 L hydrogen exposure. See also at the 'Vibrational 
properties measured by HREELS' section in 'Results'. 



Methods 

Details about SiC samples and hydrogen exposures. We use high-quality sin- 
gle-domain single- crystal 3C-SiC(001) few [im thin films grown on a vicinal (4°off) 
Si(OOl) substrate by NovaSiC and CEA-LETI. Low (10^^ atoms cm^^) and 
medium (10^^ atoms cm~ ^) n-doped samples were used for HREELS and SR-PES, 
respectively. Note that, with these two SR-PES and HREELS techniques, the 
information is collected in a 7 x 10 '^'^ jim^ area. HREELS and SR-PES experiments 
are performed at a base pressure lower than 2 x 10~ ^"^mbar. In situ atomic H 
exposures are performed with high-purity research-grade molecular H2 introduced 
by a capillary tube leading to the SiC sample with its surface maintained at 300 °C 
and dissociated by a heated tungsten filament. To avoid contamination, thorough 
baking of the gas lines and outgassing of the tungsten filaments are performed. The 
absolute determination of H coverage for experiments performed in different 
laboratories is not straight forward since it depends critically on the efficiency of 
the H2 cracker. Then, the H atoms react with the 3x2 surface reconstruction of 
the 3C-SiC(001)-3 x 2 surface (which is inert to molecular H2). However, after 
dissociation, some of the H atoms recombine to re-form H2. This recombination 
process depends very much on the geometry of the vacuum chambers, which are 
very different in the HREELS and the photoemission experimental systems. So, 
contrary to other adsorbed species, the problem with atomic H experiments is that, 
the exposure is known, not the exact coverage. Here, the 36L H PES (Synchrotron 
SOLEIL Tempo beamline) and the SOL H HREELS (IMEM-CNR & University of 
Genoa joint laboratory) correspond both to the H-coverage leading to the metal- 
lization, that is, to the calculated 8 H metallic structure, so they are comparable. All 
other details, in particular about sample preparation and atomic H/D exposures are 
available elsewhere^-i0'23,25,26,30_ 



HREELS experiments. The HREELS measurements were carried out in Genoa 
using a Delta 0.5 SPECS spectrometer in-specular and out-of-specular, at a scat- 
tering angle 9 = 60° from the surface normal, with primary electron beam energies 
£0 = 2.5 and 10 eV at 4meV energy resolution. The accuracy in vibrational fre- 
quencies determination is, however, better and estimated at ± 1 meV, that is, 
± 8 cm~ ^ The data are collected at room and low temperatures. The achievement 
of the metallic state is evidenced in the HREELS data through the elastic peak 
broadening as already shown above in Fig. 3 and also in Fig. 8 for the clean and the 
501 hydrogen exposed 3C-SiC(001)-3 x 2 surfaces, respectively, at 30, 105 and 
300 K. We notice that the specular peak broadens with temperature in both cases, 
and that the effect is much stronger for the H-covered case. This phenomenon has 
been investigated by Persson and Demuth and published in a seminal paper in 
198438 particular, they solved analytically the cases corresponding to: (a) a 
doped semiconductor: Applying the formula given in ref. 38 to the bare surface, we 
obtain a carrier density of 1,017 carriers per cm^ which compares reasonably well 
with the nominal doping of our sample of ^1,016 carriers per cm^; and to (b) a 
metallic film on a semiconducting surface. This case was solved for the limits in 
which, the plasmon energy is much smaller and much larger than thermal energy 
at room temperature. These limits are verified for the cases of Si(lll)-7 x 7 
(plasmon energy fiofp = 0.2 meV) and for ultrathin films of Au and Pd on Si(lll) 
{hcop ^ kT). In the first case, the broadening is due to the multiple excitation of the 
plasmon, whereas in the second case, it is associated to generation of electron hole 
pairs (Drude damping). For Si(lll)7 x 7 the carrier density is 1/49 per substrate 
unit cell and the effective mass of the electron 60 free electron masses. For 8 H-3C- 
SiC(001)-3 X 2, the carrier density is ~ 15 times larger. Assuming a unitary 
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effective mass, we get an upper limit for the plasmon energy of 6 meV, which is in 
between the limits evaluated by Person and Demuth^^. A quantitative evaluation of 
the broadening is therefore not possible, but one can safely derive that the observed 
temperature-dependent broadening indicates a metallic behaviour. 

Synchrotron radiation-based PES experiments. Valence-band SR-PES experi- 
ments are performed at 300 K, using the ultrahigh vacuum photoemission 
experimental station of the TEMPO beam line connected through an ondulator to 
the light emitted by the third generation synchrotron radiation SOLEIL light 
source. The diameter of the SR beam and the size of the photoelectrons collected 
area are 300 |im. Photoemission spectra are measured using a Scienta SES2002 
hemispherical electron energy analyser equipped with a fast two-dimensional delay 
line detector. The overall energy resolution combining the contributions of the soft 
X-ray beam monochromator and analyser energy selection is, in our experimental 
conditions, better than 50 meV. 

The experimental station allows direct in situ atomic H exposures during 
photoelectron spectroscopy data acquisition up to partial pressures of 10~^mbar. 
Particular care was taken in protecting the sample surface from metallic 
contamination of the surface induced by the high temperature tungsten filament 
used to dissociate the H2 molecules. 

The achievement of the metallic state is evidenced through the real-time 
measurement of the Fermi level kinetic energy region and by monitoring the 
corresponding build-up of the DOS as presented in Fig. 4. 

Ab-initio simulations. The interaction of H atoms with this surface is investigated 
also by ab-initio calculations on a model, which consists of five Si layers and four C 
layers. The top layers represent the (3 x 2) reconstructed surface. For convenience, 
the bonds of the carbon layer at the bottom of the model are saturated with 
selenium atoms. The calculations are based on DFT^^'^^ with the generalized 
gradient approximation in the form proposed by Perdew, Burke and Ernzerhof^^. 
The Kohn-Sham equations are solved with the all-electron frozen-core projector 
augmented wave method^^ as implemented in the Vienna ab-initio Simulation 
Package^^"^^ and integrated in the MedeA® computational environment^^. Locally 
stable equilibrium structures are determined by minimizing energy and forces 
using a conjugate gradient method. The topmost six layers are relaxed so that the 
forces on the atoms are <0.02eVA~ ^ The three layers at the bottom are kept 
frozen. The vibrational frequencies are computed using a finite difference method^^ 
as also implemented in the MedeA* environment^^. In these phonon calculations, 
the bottom layers are kept frozen and phonon dispersions originating from the 
vibrations of all other atoms are computed. For the clean 3C-SiC(001)-3 x 2 
surface, comparison between computed DOS and measured one by valence band 
photoemission is shown in Fig. 8, with the valence band photoemission spectrum 
taken from ref. 23. Note the remarkable agreement between theory and experiment, 
in particular for the shape and energy position of the bulk (B) and surface (S) 
electronic states. 
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